A FAMILY OF C° FINITE ELEMENTS FOR 
KIRCHHOFF PLATES I: ERROR ANALYSIS 



L. BEIRAO DA VEIGA*, J. NIIRANENt, AND R. STENBERG* 

Abstract. A new finite element formulation for the Kirchhoff plate model is presented. The 
method is a displacement formulation with the defiection and the rotation vector as unknowns and 
it is based on ideas stemming from a stabilized method for the Reissner-Mindlin model [11] and a 
method to treat a free boundary [5]. Optimal a-priori and a-posteriori error estimates are derived. 
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1. Introduction. A conforming finite element method for the Kirchhoff plate 
bending problem requires a -continuity and hence leads to methods that are rarely 
used in practice. Instead, either a nonconforming method is used or then the model is 
abandoned in favor of the Reissner-Mindlin model. For the latter, there exist several 
families of methods that have rigorously been shown to be free from locking and 
optimally convergent. 

A natural idea is to consider the Kirchhoff model as the Hmit of the Reissner- 
Mindlin model when the plate thickness approaches zero and to use a good Reissner- 
Mindlin element with the thickness (after a scaHng, see below) representing the pa- 
rameter penalizing the Kirchhoff constraint. In this approach, there are two obstacles. 
First, for a free boundary, this leads to a method which is not consistent. In the lit- 
erature, this point is often ignored since mostly the clamped case is considered. A 
remedy to this was developed by Destuynder and Nevers who showed, that the con- 
sistency is obtained by adding a term penaHzing the tangential Kirchhoff condition 
along the free boundary [5]. Even if this modification has been done, there remains a 
second drawback. In order that the solution to the penalized formulation is close to 
the exact solution, the penalty parameter should be large. This, however, leads to an 
ill-conditioned discrete system. 

Our aim in the present paper is to present a family of Kirchhoff plate bending 
elements for which the convergence rate is optimal even in the presence of free bound- 
aries. The method is a formulation combining the ideas from the stabilized method for 
the Reissner-Mindlin plates presented in [11] and the treatment of the free boundary 
presented in [5] (cf. [1, 2] as well). The family includes "simple low-order" elements 
and it is well-conditioned. In the second part [3] of this paper, we give the results of 
numerical tests. 

The paper is organized as follows. In the next section, we describe the plate 
bending problem, and in Section 3, we introduce the new family of finite elements. In 
Section 4, an a-priori error analysis is derived. This analysis leads to optimal results, 
both with respect to the regularity of the solution and to the polynomial degree used. 
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In Section 5, an a-posteriori error analysis is performed. We derive a local error 
indicator which is shown to be both reliable and efficient. 

2. The KirchhofF plate bending problem. We consider the problem of bend- 
ing of an isotropic linearly elastic plate and assume that the undeformed plate mid- 
surface is described by a given convex polygonal domain f2 C K^. The plate is 
considered to be clamped on the part Tq of its boundary dil, simply supported on 
the part Fs C dfl and free on Fp C dfl. The deflection and transversal load are 
denoted by w and g, respectively. 

In the sequel, wc indicate with V the set of all corner points in Fp. Moreover, 
n and s represent the unit outward normal and the unit counterclockwise tangent 
to the boundary. Finally, for points x € V, we introduce the following notation. We 
indicate with ni and si the unit vectors corresponding, respectively, to n and s on 
one of the two edges forming the boundary angle at x; with n2 and S2 we indicate the 
ones corresponding to the other edge. Note that which of the two edges correspond 
to the subscript 1 or 2 is not relevant. 

The classical Kirchhoff plate bending model is then given by the biharmonic 
partial differential equation 



DA'^w = g inO, (2.1) 



the boundary conditions 



«' = 0, 1^=0 onFc, 

w = 0, n- Mn = on Fs , (2.2) 

n ■ Mn = , ^ (s • Mn) + (div M) ■ n = on Fp , 
and the corner conditions 

(si • Mni){x) = (s2 • Mn2){x) Va; e V. (2.3) 

Here 

is the bending rigidity, with E, being the Young modulus and the Poisson ratio for 
the material, respectively. Note that for the shear modulus G it holds 

The moment tensor is given by 

M{Vw) = D((l - i/)e(Vw) + zvdiv(Vw) J) , (2.6) 

with the symmetric gradient e, and the shear force by 

Q = -divM. (2.7) 

Note, that the independence of the Poisson ratio v in the differential equation (2.1) is 
a consequence of cancellations when substituting (2.6) and (2.7) into the equilibrium 
equation 

-diyQ = g. (2.8) 
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For the analysis below, it will be convenient to perform a scaling of the problem by 
assuming that the load is given by g = Gt^f, with / fixed. Then the differential 
equation (2.1) becomes independent of the plate thickness: 

^ A'^w = f, infl. (2.9) 



6{l-u) 

Furthermore, we use the following scaled moment tensor m given by 

M (Vw) = Gi^m(Vw) (2.10) 
and the shear force q is defined by 

Q = Gt^q. (2.11) 

The unknowns in our finite element method will be the approximations to the 
deflection and its gradient, the rotation f3 = 'Vw. With this as a new unknown, our 
problem can be written as the system of partial differential equations 

Vw-/3 = 0, (2.12) 

-divq = /, (2.13) 

Lf3 + q = 0, inO, (2.14) 

the boundary conditions 

w = 0, /3 = 0, on Tc , (2.15) 

w = 0, /3-s = 0, n- m{/3)n = , on Fs , (2.16) 

— - /3 • s = , n- m{(3)n = , ^ (« ' m{l3)n) - q • n = , on Fp , (2.17) 
and the corner conditions 

{si-m{f3)ni){x) = {s2-m{f3)n2){x) Va;eV. (2.18) 
The operator L is deflned as 

I,/3 = divm(/3) , (2.19) 
and the scaled bending moment is considered as a function of the rotation: 

m(/3) = ^(£(^) + ^div/3J). (2.20) 

In the sequel, we will often write m instead of m{/3). We further denote 

a(/3,r7) = (m(/3),£(r/)). (2.21) 

In order to neglect plate rigid movements and the related technicalities, we will 
in the sequel assume that the one-dimensional measure of Fc is positive. 
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3. The finite element formulation. In this section, we will introduce our finite 
element method. Even if our method is stable for all choices of finite element spaces, 
we will, for simplicity, present it for triangular elements and for the polynomial degrees 
that yield an optimal convergence rate. Hence, let a regular family of triangular 
meshes on O be given. For the integer fc > 1, we then define the discrete spaces 

Wh = {w&W \ w\K& Pk+i{K) VK e Ch} , (3.1) 

Vh = {v€V I T7|K G [Pk{K)]^ yK e Ch} , (3.2) 

with 

W={v€H^{n)\v = OonTcUTs} , (3.3) 

V ={r]G [H^{il)f I Tj = on Tc , T7 • s = on Ts} . (3.4) 

Here Ch represents the set of all triangles K of the mesh and Pk (K) is the space of 
polynomials of degree k on K. In the sequel, we will indicate with hx the diameter 
of each element K, while h will indicate the maximum size of all the elements in the 
mesh. Furthermore, we will indicate with E a general edge of the triangulation and 
with He the length of E. The set of all edges lying on the free boundary Fp we denote 
by J^h- 

Before introducing the method, we state the following result which trivially follows 
from classical scaling arguments and the coercivity of the form a. 
Lemma 3.1. There exist positive constants Cj and Cj such that 

Ci ^ hULci>\\lK<a{<P,<P) y<P&Vh, (3.5) 
KeCh 

C'j Y,hE \\mnsm\lE < a{(p, ct>) V0 e ^ , (3.6) 

where the operator mns{4>) = s ■m{4))n with n, s, being the unit outward normal and 
the unit counterclockwise tangent to the edge E, and with m defined in (2.20). 

Let two real numbers 7 and a be assigned, 7 > 2/Cj and < a < Ci/4. Then, 
the discrete problem reads as follows. 

Method 3.1. Find (wh,f3,J e Wh x- Vh, such that 

Ah {wh,^h- V, V) = if, v) ^{v, v)^WhxVh, (3.7) 
where the form Ah is defined as 

•Ah{z, 0; V, rj) = Bh{z, 4>; v, rj) + Vh{z, 0; v, rj) , (3.8) 

with 

Bh{z, (j); V, rj) = a{(l), ry) - ^ ahj^{L^, Lr))K 

+ ^ —-^{Vz-(f)-ahjiL4>,Vv-r]-ahjiLr])K (3.9) 

and 

T^h{z, 0; V, T]) = (m„s(0), [Vv - rj] ■ s)rp + {[Vz - <p] ■ s, m„s(»7))rF 

+ H TT-^t'^^ - 0] • [V^; - V] ■ s)e (3.10) 
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for all {z , (p) , {v , ij) € Wh x V?,. Here (•, •)rF '^"'^ {'j')e denote the -inner products 
on T-p and E, respectively. 

The bilinear form Bh constitutes the Reissner-MindHn method of [11] with the 
thickness t set equal to zero, while the additional form T>h is introduced in order to 
avoid the convergence deterioration in the presence of free boundaries. 

Furthermore, we introduce the discrete shear force 

Qh\K = ^("^"^h -f3h- ahlLI3f,)\K yK&Ch. (3.11) 
We note that, due to (2.14) and (2.12), it holds 

q\K = ^(^^ - /3 - Oih%L(3)\K "^K&Ch, (3.12) 

and hence it follows that the definition (3.11) is consistent with the exact shear force. 

For simplicity, in the rest of this section we assume that the defiection w belongs to 
H^{Q); this is a Yevy reasonable assumption, as discussed at the end of this section. 
Note as well that, with some additional technical work involving the appropriate 
Sobolev spaces and their duals, such assumption could be probably avoided. The 
following result states the consistency of the method. 

Theorem 3.2. The solution {w,(3) of the problem (2.14)-(2.18) satisfies 

Ah{w,f3;v,'n) = {f,v) y{v,'n) GWhxVh. (3.13) 

Proof. The definition of the bilinear forms in Method 3.1, recalling (2.14) and the 
expression (3.12), give 

Bh{w,f3;v,r]) = a{f3,r]) - ^ ah]i{L(3,Lr^)K 

+ ^ ^2 — P — ah%Ll3, Vv — rj — ah\L'q)K 

= a(/3,r7)+ ^ ah\{q,Lr})K + ^ {q,Vv - r] - ah\Lr})K 

KeCh KeCh 
= a{/3,rj) + {q,Vv-rj). (3.14) 

First, by the definition (2.21), then integrating by parts on each triangle, finally using 
the regularity of the functions involved, and the boundary conditions (2.15), (2.16) 

on Fc, Fs, respectively, we get 

a{f3, V) + (q, Vw - rj) = (m(/3), £(77)) + (q, Vw - q) 

= -{L(3 + q,ri) + {m{f3) ■ n, T])rp - (div q, v) + {q ■ n, v)rp ■ (3.15) 

Recalling (2.14) and (2.13), the identity above becomes 

ail3, r?) + {q, Vv - rj) = (/, v) + {m{fi) ■ n, rt)r^ + {q-n, v)r^ , (3.16) 

while, using the boundary conditions of (2.17) on Fp and integration by parts along 
the boundary, finally leads to 

a(/3, rj) + {q, Vv-rj) = (/, v) - (m„«(/3), [Vv - r,] ■ s)rp • (3.17) 
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Due to (2.17), we have 

'Dh{w,l3;v,r]) = (m„s(/3), [Vw - t/] • s)rp + ([Vw - /3] • s,mns{'n))rF 
+ T^ii^w - /3] ■ s, [Vu - rj] ■ s)e 

= (m„,(/3),[Vw-r7] •s)rp . (3.18) 

The result now directly follows from (3.14), (3.17) and (3.18). □ 

Remark 3.1. Ij the Reissner-Mindlin method of [11] without the additional form 
Vh is employed by setting t = 0, then in the presence of a free boundary we obtain 

Bh{w,/3;v,ri) = {f,v) + {mns{l3),[^y-v]-s)r^ V(t;, rj) e x ^ . (3.19) 

Therefore, this would lead to an inconsistent method. We return to this in Remark 
4-1 below. 

4. Stability and a-priori error estimates. For {v,ri) G W/(X V/j, we introduce 
the following mesh dependent norms: 

\{v,v)\l= ''kW'^^-vWIk, (4.1) 
KeCh 

MIh = \\v\\1+ E \-\Ik+ E ^^^'1111^1110.^;+ E ''e'W^WIe, (4.2) 

lll(^,»7)IIU = l|r7l|i + lklk. + |(^;,»7)|/., (4.3) 

where |.| represents the jump operator and Xh denotes the edges lying in the interior 
of the domain fl. 

In [9], the following lemma is proved. 

Lemma 4.1. There exists a positive constant C such that 

\\v\\2,h<C{\\r]\\i + \\v\\i + \{v,r))\h) V(t;, r?) G x . (4.4) 

Using the Poincare inequality and the previous lemma, the following equivalence 
easily follows: 

Lemma 4.2. There exists a positive constant C such that 

C|||(^;,r/)||U< ||»7||i + |(t;,»7)U< |||(^^,r?)||U y{v,v) & x Vh . (4.5) 

We now have the following stability estimate. 

Theorem 4.3. Let < a < Cj/A and 7 > 2/Cj. Then there exists a positive 
constant C such that 

Ah{v,ii;v,n)>C\\\{v,ii)\\\l \l{v,ii)€WhxVh. (4.6) 
Proof. Using the first inverse estimate of Lemma 3.1 we get 

> (i-^)«(^>^)+ E -rrllV^^-^-"/^Ki'?llo,if (4.7) 



A family of C" finite elements for Kirchhoff plates I: Error analysis 



7 



Next, using locally the arithmetic-geometric mean inequality with the constant 7//1B, 
then the second inverse inequality of Lemma 3.1, we get 

Vh{v,'n;v,r)) 

= J2 {^mns{r)),[^v-il]-s)E + ^\\[yv-ri]-s\\lE) 

= - E 7"^^B l|TOrxs(»7)llo,i<; 

7-^ 1 

>-^^a{r],r,)>--a(r],T,). (4.8) 

Joining (4.7) with (4.8) and using Korn's inequality we then obtain 
^hiy, V, V, ri) + Vh{v, T]; v, rj) 

> C{\\r)\\l + ^ ^llVt^ -ri- ahlLriWl^) . (4.9) 
KeCh ^ 

From the triangle inequality, again the inverse estimate of Lemma 3.1 and the bound- 
edness of the bilinear form a it follows 



E -^\\^v-v-ahj,Lrj\\lj,+ ^ ahl\\Lri\\lj,) 
E ^\\'^^-'n-(^hlLr)\\lj, + a{r,,ri)) 

E ^rW'^^-n-ahlLnWlj. + WnWl), (4.io) 

KeCh ^ 

which combined with (4.9) gives 

Ahiv,rj;v,rj)>C{Ml + \iv,ii)\l). (4.11) 

The result then follows from the norm equivalence of Lemma 4.2. □ 

We can now derive the error estimates for the method. We note that the assump- 
tions of the theorem are supposed to be valid for the further results below as well, 
hence not repeated in the sequel. 

Theorem 4.4. Let < a < Cj/A and 7 > 2/C|. Let {w,(3) be the exact solution 
of the problem and let {wh,/3f^) be the approximate solution obtained with Method 3.1. 
Suppose that w G if*+^(f2), with 1 < s < k. Then it holds 

\\\iw- WH,0- /3JIIU < Ch^\\w\U+2 . (4.12) 
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Proof. Step 1. Let {wi,l3j) G Wh x Vh be the usual Lagrange interpolants to 
w and f3, respectively. Using first the stability result of Theorem 4.3 and then the 
consistency result of Theorem 3.2 one has the existence of a pair 

iv,r))eWhxVh, |||(^;,r,)||U<C, (4.13) 

such that 

\\\{wh-wi,f3h- /3i)\\\h < Ahiwh-wi,f3h-(3j;v,'n) 

= Ah{w-wi,P-f3j;v,r]), (4.14) 

where we recall that Ah = Bh + T)^ . 
Step 2. For the Bh-part, we have 

Bh{w-wi,0-(3i;v,r])=ai0-f3i,r])- ^ a/ii(X(^ - ^9,), Lt?)^ (4.15) 

KeCh 

+ Yl - - ('^ - -^z) - c^I^kHI^ - /3/), Vv-Tj- ah\Lri)K • 

Due to the first inverse inequality of Lemma 3.1, we get 

( ^ hlWLnWl^y' <C\\\{v,vmh (4.16) 

and 

( E -^l|V^-»7-a/iiir,||2_^)'^'<q||(^,,7)||U. (4.17) 

Using these bounds in (4.15) and recalling (4.13) we obtain 
Bh{w -wi,(3- f3j;v,r]) 

<c(\\\{w-wj,f3-(3j)\\W + { hlW-f^AlKf'^). (4.18) 

K&Ch 

Substituting the definition of the norm (4.3) in (4.18), using the triangle inequality, 
and finally applying the classical interpolation estimates it easily follows 

Bh{w -wi,/3- Pi; v, r,) < Ch'{\\w\U2 + ||/3|U+i) . (4.19) 

Step 3. For the 'D/j-part in (4.14), we have, by the definition (3.10), 

T>h{w -wi,/3-(3j; V, rj) = {mns{0 - /3j), [Vi; - 77] • s)rF 
+{[V{w - wi) -{13- (3j)] ■ s, m„s(T7))rp 

=: T1+T2 + T3. (4.20) 
Scaling arguments give 

II - V] ■ s\\Ie < II Vt; - vWIe < Gh-j,\E) W^^ ' ^llW) (4-21) 
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for all E G J^h, where K{E) is the triangle with E as an edge. The Z^-Cauchy-Schwartz 
inequality, the bound (4.21) and the norm definition (4.3) now give 

1/2 / ^ . , _ , xl/2 



Tl < ( ^ hK(E)\\mns{P - f3l)\\lE) ( E f^K\E)\\i^^-V]-4lE) 

1 /2 

<C( E ^^(^^)ll"^-(/3-/3/)llo.i=;) \\\{v,ri)\\\h. (4.22) 



EeJ'h 

Recalling the bound (4.13), classical polynomial interpolation properties give 

1/2 

hK(E)\\mns{l3 - (il)\\o,Ej 

' EeJ^h 

Again, by scaling we have 



1/2 

T^<C(^^ hKiE)\\mns{p-f3i)\\lE) <Ch'mU+,. (4.23) 



\\mns{v)\\lE < h-^\E)\v\lKiE) ^EeTh. (4.24) 
The Z^-Cauchy-Schwartz inequality, this bound and the norm definition (4.3) give 

T2<[ E h-^\E)\\^^^-^l)-il^-(^l)\\lEY\ E hKiE)\\mnsmlE)"^ 

EeJ^h EeJ^h 

<C( E Mi^)l|V(«'-t"/)-(/3-/3,)||^,^)'^'|||(t;,»7)||k. (4.25) 
Ee:Fh 

Recalling the bound (4.13), classical polynomial interpolation estimates give 

1 /2 

r2<c( E h-^\E)\\'^i^-^l)-((^-(^l)\\lE) 

<C/i^(||/3|U+i + ||«;||.+2). (4.26) 
The bound for T3 follows combining the same techniques used for T\ and T2; we 

get 

Ti<Ch'{\\l3\\s+x + \\w\\s+2). (4.27) 
Now, joining all the bounds (4.14), (4.19), (4.20), (4.23), (4.26) and (4.27) we obtain 

\\\{wh-wi,(3^-(3Mh < Ch^iMs+i + \Ms+2) . (4.28) 

The triangle inequality and the classical polynomial interpolation estimates (recalling 
that /3 = Vw) then yield 

\\\{w - Wh, f3 - f3h)\\\H < Ch' {mU+i + \Ms+2) <Ch'Ms+2. (4.29) 

Note that the result holds for real values of the regularity parameter s since the 
interpolation results used above are valid for real values of s. □ 

Remark 4.1. As noted in Remark 3.1, the limiting Reissner-Mindlin method 
(i.e., without the additional correction Vh) is inconsistent. Regardless of the solution 
regularity and the polynomial degree k, the inconsistency term can only be bounded with 
the order 0{h^^^). As well known (see for example [10]), the inconsistency error is a 
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lower bound for the error of finite element methods. As a consequence, the numerical 
scheme will not converge with a rate better than h^^^ «/ Tp 7^ 0- This observation is 
also confirmed by the numerical tests shown in [3]. See [4] for other numerical tests 
regarding this issue. Note further, that this boundary inconsistency term is connected 
not only to the formulation in [11] but is common to any other Kirchhoff method 
which follows a "Reissner-Mindlin limit " approach. 

For the shear force, the practical norm to use is the discrete negative norm 

||r||_i,, = ( ^ hl\\r\\l^f\ (4.30) 
KeCu 

Since we assume that w G with s > 1, we have q G [L^(n)]^ and from the 

estimates above it immediately follows: 
Lemma 4.5. It holds 

\\q-qd-i,h<Ch'\\w\U2. (4.31) 



From this it follows a norm estimate in the dual to the space 

= {r/ G [H\n)]'^ I 77 = on Fc, r? • s = on Fp U Fs} , (4.32) 

i.e., in the norm 

||r||_i,* = sup (4.33) 
riev. Mil 



The result we have is the following. 
Lemma 4.6. It holds 



\\q-qh\\-i,* <Ch'\\w\\s+2. (4.34) 



Proof. The proof is essentially an application of the "Pitkaranta-Vcrfiirth trick" 
(see [8, 12]). By the definition of the norm || • there exists a function 77 G V* 

such that 

||g-'7,||-i,*<('7-'7„»7), \\ri\\i<C. (4.35) 

Using a Clement type interpolant we can find a piecewise hnear function rjj G V* 
such that it holds 

h'K^\\r]-rij\\s,K <C\\'n\\i,K <C' , s = 0, 1 (4.36) 

for all K G Ch. Using the Cauchy-Schwartz inequality, the bound (4.36) with s = 
and the definition (4.30) it follows 

(g - Qh^ 'n) = {q- q^, r)-'qj) + {q- qf,, Vi) 

<C\\q-qf,\\.i,H + {q-qH,Vi)- (4.37) 

Note that r/j is both in Vh and V* ; moreover Lrjj = on each element K of Ch. As 
a consequence, using (3.7), (3.11), (3.12) and Theorem 3.2, it follows 

(9 - Qh^ Vi) = - f^h, Vi) + {[^Wh - Ph)] ■ s> Afns(»7/)>rF 

=: T1+T2. (4.38) 
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Due to the continuity of the bihnear form and using bound (4.36) with s = 1 it 
immediately follows 

Ti < C\\I3 - /3J|i < C\\\{w -Wh,l3- p^)\\\h . (4.39) 

Using first the Cauchy-Schwartz inequality, then the Agmon inequality, finally the 
bound (4.36) with s = 1, Lemma 3.1 and the definition (4.3), we get 

T2< ( ^ hl'\\Vwh-l3^)\\lEY\ E hE\\Mns{rij)\\lj,)"^ 

K&Ch 

<C\\\{w-Wh,l3-l3h)\\\h, (4.40) 

where in the last inequality wc implicitly used the relation — /3 = 0. Combining 
(4.35), (4.37) with (4.38), (4.39) and (4.40), it follows that 

\\q - QhW-i,* < C{\\q - qA-i,h + -Wh,/3- /3J||U) . (4.41) 

Joining (4.41), (4.31), and using Theorem 4.4 the proposition immediately follows. □ 
The regularity of the solution to the Kirchhoff plate problems for convex polygonal 

domains, with all the three main types of boundary conditions, is very case dependent. 

We refer for example to the work [7] in which a rather complete study is accomplished. 

Note that, if / G H~^{^1), in most cases of interest, the regularity condition w e 

H^{iY) is indeed achieved. 

Note further, that with classical duality arguments and technical calculations it 

is possible to derive the error bound 

\\w-Wh\\i<Ch'+^\\w\U2, (4.42) 

if the regularity estimate 

ll«^l|3<C||/||_i (4.43) 
holds. Moreover, if fc > 2 and the regularity estimate 

\HU < C\\fh (4.44) 

is satisfied, then it holds 

\\w - WhWo < Ch'^^\\w\\s+2 . (4.45) 

5. A-posteriori error estimates. In this section, we prove the reliability and 
the efficiency for an a-posteriori error estimator for our method. To this end, we 
introduce 

:= h%\\f + divq,,\\li, + h],^\\Vw,, - pj^j, , (5.1) 
Vl ■■=h%\\lq^-nl\\lE + hE\\lm{Pf,)nl\\lE, (5.2) 

V^E ■= ^i5|l™mi(/3ft)llo,_E : (5.3) 

d 

Vf,E ■■= hE\\mnn{fih)\\o,E + ^lll^™ns(^ft) " 9?, " "-llo.B > (5-4) 
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where He denotes the length of the edge E and [•] represents the jump operator 
(which is assumed to be equal to the function value on boundary edges). Further, for 
a triangle K ^ Ch we denote the sets of edges lying in the interior of on Fs and on 
Fp, by liK), S{K) and F{K), respectively. By Sh we denote the set of all edges on 
Fs and by Xh the ones lying in the interior of the domain. 

Given any element K &Ch, let the local error indicator be 

Vk '■={flK + \ J2 Yl ^S,E + Yl "^f^e) ' > (5-5) 

EeI{K) EeS{K) EeF{K) 

Finally, the global error indicator is defined as 

1 /2 

^ ■= ( E • (5-6) 

KeCh 

Remark 5.1. It is worth noting that, by the definition (3.11), 

{q^ + Lp^)\K=J^{^Wh-l3f,)iK yKeCh, (5.7) 

which is the reason why there appears no terms of the kind \\qf^ + Lf3f^\\o,K in the error 
estimator. We note as well that scaling arguments give 

Yl hE'\\Vwh-f3X,E<C Y hK^W'^^h-Mlx, (5-8) 

which is the reason why there appears no boundary terms of the kind WVwh— f3h\\o,E- 

5.1. Upper bound. In order to derive the reliability of the method we need the 
following saturation assumption. 

Assumption 5.1. Given a meshCh, letCh/2 be the mesh obtained by splitting each 

triangle K E Ch into four triangles connecting the edge midpoints. Let {wh/2- 

be the discrete solution corresponding to the mesh Ch/2. We assume that there exists 

a constant p, < p < 1, such that 

\\\{W - Wh/2,f3- f3h/2)\\\h/2 + \\q - 9/1/211-1,* 

< p{\\\{w -Wh,f3- f3^)\\\h + Wq-QhW-i.*) , (5-9) 
where by ||| • \\\h/2 we indicate the mesh dependent norm with respect to the new mesh 

Ch/2- 

In the sequel, we will need the following result. 
Lemma 5.1. Let, for v G Wfi/2, the local seminorm be 

\v\2,h/2,K={ Y MIk')'^^- (5-10) 

Then, there is a positive constant C such that for all v € there exists vj € Wh 

with the bound 

\\v - viWo^K + H^k'^Wv - viWo^dK <ChWv\2^h/2,K yKeCh- (5.11) 
Moreover, vj interpolates v at all the vertices of the triangulation Cy^i2- 
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Proof. We choose Vj as the only function in H^{Vl) such that 

vi\K e P2{K) yKeCh, 

vi{x) = v{x) yx e Vh/2 , (5.12) 

where Vh/2 represents the set of all the vertices of 0^/2- Note that it is trivial to check 
that Vi G Wh for all fc > 1. Observing that 

\v\2,h/2,K + vGWh/2,KGCh, (5.13) 

is indeed a norm on the finite dimensional space of the functions v e Wh/2 restricted 
to K, the result follows applying the classical scahng argument. □ 

For simplicity, in the sequel we will treat the case Fs = 0, the general case 
following with identical arguments as the ones that follow. We have the following 
preliminary result. 

Theorem 5.2. It holds 

\\\iwh/2-Wh,f3h/2-f3h)\\\h/2<Cr,. (5.14) 

Proof. Step 1. Due to the stability of the discrete formulation, proved in Theorem 
4.3, there exists a couple {v,r}) G Wh/2 x Vh/2 such that 

|||(t^,r/)||U/2<C (5.15) 

and 

\\\{Wh/2-Wh,l3h/2- l^h)\\\h/2 < Ah/2iWh/2 - Wh, f3h/2 - 0h^v,r]) . (5.16) 
Furthermore, we have 

Ah/2{wh/2,f3h/2'-:'",V) = {.f,v) ■ (5.17) 

Step 2. Simple calculations and the definition (3.11) give 

Bh/2{wh,0h; V, rj) = a{l3h, r]) - ^1 ^if^KiL^h, Lr])K 

+ X] - Ph- Oih\Ll3h , - 77 - ahjiLr])K 

= a{f^h^v)- i^Wh - l3h,Lr])K + ^ (qh,^v-r])K 

KeCh/2 KeCh/2 
+Ri{wh,(3h;v,r]) 

= Bh{wh,Ph, V, ri) + Ri {wh, /3h; v, rj) , (5.18) 
where Qh is defined as in (3.11), i.e., based on the coarser mesh, and 

Ri{wh,^h-iV,r)) = Y -rri'^'^h- /3h,'^v -r])K 

KeCu/2 



Y -^V^Wh-^h,^v-v)K. (5.19) 
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The last term on the right hand side is well defined since Vv—r] is piecewise i^-regular. 

Let now indicate the set of all edges of 0^/2 lying on Fp. Adding and 

subtracting the difference between the two forms it then follows 

1^/1/2 (wh, /3^; V, vi) = Vh{wh,l3h, V, rj) + R2{wh,f3h; v, rj) , (5.20) 

where 

-R2(wh,/3^;v,r7) = V -^{[Vwh - (3,,] ■ s,[Vv - rj] ■ s)e 

- E ^{[^Wh-l3h]-s,[^v-ri\-s)E, (5.21) 

and where the first member on the right hand side is indeed well defined due to the 
piecewise regularity of {v,ri). We will denote 

R{wh,0h;v,'n) = Ri{wh,f3h;'v,'n) + R2{wh,fih'V,'n) . (5.22) 

Joining (5.17)-(5.21) then yields 

•AhMwh,f3hi ^> = -Ahiwh^Ph^ V, r}) + R{wh,l3h; v, rj) . (5.23) 

Step 3. Let vj G Wh be the interpolant defined in Lemma 5.1 and let rjj € Vh be 
the piecewise linear interpolant to rj. First, we have 

Ah{wh,f3h; vi,r]i) = (/, vi) . (5.24) 

This, together with (5.17) and (5.23) gives 

A/2(wh/2 - Wh, /3/1/2 - (^h^V, Tj) 

= ■Ah/2{u!h/2,f3h/2;i>,v) - ■Ah/2(wh,l3h;v,r]) 

= ■^h/2{wh/2,f3h/2l'",v) - ■Ah{wh,(ih'^v,r]) - R{wh,Ph'^v,r]) 

= {f,v- vi) - Ah{wh,Ph^ v-vi,r)- rjj) - R{wh,/3h; v, 77) . (5.25) 

Step 4- Next, we bound the last terms above. Recalling that Cfi/2 is a subdivision 
of C/i, the Cauchy-Schwartz inequality, (4.3) and (5.15) give 

\Ri{wh,f3h;v,'n)\<2\ ^ _l_(v«;;,-/3^,Vf-»7)x| 
KeCh/2 " ^ 



E ^W^^^-MIk)'^' ■ (5.26) 



KeCh/2 



Using scaling and arguments similar to those already adopted in (5.26) it can be 
checked that 

\R2{wH,0h;v,v)\<C{ Yl ^W'^^h - MIk)'^' ■ (5-27) 
Kec,,2 '^^ 



A family of C° finite elements for Kirchhoff plates I: Error analysis 



15 



Combining (5.26) and (5.27) we get 

\R{wh,/3h;v,v)\ < \Riiwh,l3h;v,ri)\ + \R2{wh,l3h;v,ri)\<Cr]. (5.28) 



Step 5. Next, we expand, substitute the expression (3.11) for Qf^ and regroup the 
terms: 



{f,v- vi) - Ah{wh,f3h;v - VI, T]- r]j) 

= {f,v -vi) - [a{l3h,Tl -Vi) - Yl (^^KiLf^h^Hv -Vi))k 



KeCh 



+ 5Z ^^(^^'^ - l^h- ah%L^f^,V{v - vi) - (tj - tjj) - ah%L{r) - rij))j^ 

+(m„s(/3ft), [V(w - vi) - (r? - jy^)] • s)^^ + {[Vwh - (3^] ■ s, rrinsiv - '7/))rp 
+ E T-d'^'^h - Ph] ■ Niv - vi) - (r? - Hj)] ■ s)^ 

KeCh 

+ (gh,V(i;-w/) - iv-vi)) 

+{mns{f3h), [V(t; - vi) - (r) - t/j)] • s}^^ + {[Vwh - Ph] ' rnns{r} - T7/))rp 



.IK 



- E T^([Vt«/.-;9ft]-5,[V(^;-t;/)]-5>^} 

-{mnsiPh),[n-r]i] ■ s)rp + ([Vw?s - /3/,] • s, m„s(r? - T?/))rp 

- E ^([Vt"/.-/3/.]-5,[r,-r7,]-s)^} 

=:A-B. (5.29) 



5iep 6. In the part A above, integration by parts, and using the fact that v{x) = 
Vi{x) at the corner points x €V, yields 



(/, v-vi)- (q^, W{v - vi)) - (mnsiPh), l^i^ - vi)] ■ s) 



d 

= (/ + div q^, v-vi) + (g^ninsiPh) -Qh '^^'"- ^J')rp • (5-30) 



The separate terms are then estimated as follows, using the Cauchy-Schwartz inequal- 
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ity and Lemma 5.1, 

\{f + diyqh,v -vi)\ = \ ^ {fh + dW qh,v - Vi)k\ 

KeCh 

< ( E hUf + diyqXK^'i E ^^K\\v-Vi\\lKy' 
KeCh KeCh 

1/2. ^ , .1/2 



< c 



E hUf + diyqJo,K) ( E I 
KeCu KeCh 

1/2 



and 



<C(E^^) (5-31) 



d d 
\{-Q^rnnsif3h) - Qh ■ n,v - vi)^J = \ E ig^^^^'^f^h) - Qh ■ n,v - vi)e\ 

f) \ 1/2 / \ 1/2 



^ ( E - 9/. • w|io,ij) ( E '^B^ii^ - ^^IIo.b) 

<^(E^M'''(E HW)' 



<C( E (5.32) 
The last term in A is readily estimated by scaling estimates and Lemma 5.1 

1/2/ ^ , , \l/2 



<C( E ^^'I|V«'/^-/3JIo,k)'^'(E ^^'ll^-^/llo,i<;)'^' 

<C^( E '^^)'''( E i-iW)'^' E '^^)'^'- (5-33) 

KeCh KeCh KeCh 

Collecting (5.30)-(5.1) we obtain 

\A\<Cr]. (5.34) 

Step 7. We will now estimate the term B. The following terms are directly 
estimated as the similar terms above 

{[S7wh-(3h]-s,mn.iv-rii)}^ +1 E ■^{i^'J^h - Ph] ' ^dv - vA ■ s) 

' Ee^h 

<Cr]. (5.35) 
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Since rii is piecewise linear, it holds Lrjj^i^ = 0. The inverse estimate then gives 

KeCh KeCh 
KeCh 

< Ct] , (5.36) 

where we in the last step used (5.7). The final step in estimating the term B is to 
integrate by parts, use the Cauchy-Schwartz inequality, interpolation estimates and 
again (5.7): 

= I - X! (LPh + Qh^V-Vi) + J2 il'^(f^h)nlv-rii)E 
KeCh Eeih 

+ {m„n{f3h), [v - Vi] ■ '^)rsurp | 

< Yl \\I^l3h + Qh\\o,K\\v-Vi\\o,K+ Yl Wl-mif^hMWoMlv-VihE 
KeCh Eeih 

+ Y N"n(/3/»)llo,£||»7-^/l|o,B 

EeShUJ^h 

< Cri . (5.37) 
Collecting (5.35)-(5.37) we obtain 

\B\<C'q. (5.38) 

Step 8. The asserted estimate now follows from (5.16), (5.25), (5.28), (5.29), 
(5.34) and (5.38). □ 

We also have the following lemma for the shear force: 
Lemma 5.3. It holds 

hh/2 - QhW-i,* < C( II IK/2 - Wh,Ph/2-l3h)\\\h/2+v) ■ (5.39) 



Proof. We start by observing that, referring to the definition (3.11) and its "/i/2" 
counterpart, Qf^ and are defined on different meshes and therefore with different 

coefRcients. However, recalling that the size ratio between the two meshes is 
bounded, it is easy to check that an opportune splitting and the triangle inequality 
give 

\\Qh/2-'lh\\-l,h<C(^ J2 \\^iWh/2-Wh)-{f3^/2-0h)\\lK 

KeCh/2 

+ Y,\\^Wh-MlK+ E hl\\L(3^/2-LMlj,). (5.40) 
KeCh KeCh,2 

The first and the last term in (5.40) can be bounded in terms of the ||| • |||/i/2 norm, 
simply using the definition (4.3) and the inverse inequality 



(5.41) 
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Therefore, recaUing the definition (5.1), we get 

\\Qh/2 - Qh\\-i,h < C(|||(w;,,/2 - Wh,f3h/2 - l3h)\\\h/2 + v) ■ (5.42) 

The transition from the \\qfi/2 ~Qh\\-'i-,h norm to the \\qh/2 norm is accom- 

plished by using the "Pitkaranta-Verfiirth trick" with steps almost identical to those 
used in Lemma 4.5, therefore omitted. □ 

Joining Theorem 5.2 and Lemma 5.3 gives the following a-posteriori upper bound 
for the method. 

Theorem 5.4. It holds 

\\\{w-Wh,0-l3h)\\\h+\\q-qh\\-i,* < Cv ■ (5.43) 

Proof. Theorem 5.2 combined with Lemma 5.3 trivially gives 

\\\iwh/2 - Wh, f3h^2 - f3h)\\\h/2 + \\qh/2 - QhW-h* < (5.44) 
From the saturation assumption it follows 

\\\{w-Wh,l3-f3,,)\\\h/2 + \\q-qh\\-i,* 
< Y^p{\\\{wh/2-Wh, 13^/2 -(3h)\\\h/2 + \\qh/2 - QhW-h*) (5.45) 

and hence the assertion follows from (5.44). □ 

5.2. Lower bound. In this section, we prove the efficiency of the error estima- 
tor. Given any edge E of the triangulation, we define loe as the set of all tlio trianglos 
K &Ch that have E as an edge. Given any K £ Ch, we define ojk as the set of all the 
triangles in Ch that share an edge with K. We then have the following lemma [6]: 

Lemma 5.5. Given any edge E of the triangulation Ch, let Pk{E) he the space of 
polynomials of degree at most k on E. There exists a linear operator 



■ Pk{E) HI{ue) (5.46) 
such that for all pk € Pk{E) it holds 

CMIe < {Pk,IiEiPk))E < WPkWlE, (5.47) 

||ni5(pfe)||o,c.^ < C2h'f\\pk\WE , (5.48) 



where the positive constants Ci above depend only on k and the minimum angle of the 
triangles in Ch ■ 

Next, we define a local counterpart of the negative norm defined in (4.33) for the 
shear force. 

||r||_i,*,a,. = sup p^. (5.49) 
nev, mi 

77=0 in Q\ojK 

We then have the following reliability result: 
Theorem 5.6. It holds 

m < C{\\\{w-Wh,(3-l3h)\\\h,u.K + \W-Qh\\-h*,^K+h]i\\f-fh\\o,u.K), (5.50) 
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where fh is some approximation of the load f. Here ||| • ||U,cj/f and \\ ■ ||o,wk represent, 
respectively, the standard restrictions of the norms \\ \ ■ \\\h and \\ ■ ||o to the domain 

UJk- 

Proof. The proof of the theorem consists of bounding separately all the addenda 
of rjK in (5.5). 

Step 1. We first bound the terms of 77^ in (5.1). Considering the right hand side 
of (5.50), the triangle inequality immediately shows that it is sufScient to bound the 
term h'^nWfh + divq^||o,A:. 

Given any K ^ Ch, let bx indicate the standard third order polynomial bubble 
function on K, scaled such that ||^_R:||L°°(ii:) = 1- Given K S Ch, let now tpx G Hq{K) 
be defined as 

ipK = {fh + divqh)b% . (5.51) 

The standard scaling arguments then easily show that 

\\fh + diyqX,K < Cifh + diyq^,^K)K, (5.52) 
\Mo,K < C\\fh + diyqJo,K ■ (5.53) 

For the first term in f]'^, the equihbrium equation (2.13) and integration by parts 

give 

/i|fllA + divq^||o,K < Ch\{fh + d:Yvqh,ipK)K 
= Chl^iif + d\yq^,LpK)K + [fh - /, Vk)k) 
= C/i|((-divq + divq^,(^_R-)/f + [fh - /, 'Pk)k) 
= Chi {{q^ - q, V^k)k + {fh - f, ^k)k) ■ (5.54) 

We note, in particular, that \7ifj^ g and \7(pK = in Q\K. Therefore, the duality 
inequality and the Cauchy-Schwartz inequality, followed by the inverse inequality and 
the bound (5.53) lead to the estimate 

Chl{{q^ - q, V^k)k + {fh - /, m)) 

< C\\q - qh\\-i,*,K h1\\'^ipK\\i,K + Ch%\\f - fh\\o,K\\v>K\\o,K 

< C{\\q - qd-i,*,K + h%\\f - fh\\o,K) \\fh + divqJ|o,if • (5.55) 

Combining now (5.54) with (5.55) gives 

hUfh + div qdo,K < C{\\q - qJ-i,.,K + hl\\f - M|o,x) • (5.56) 

The second term of f]1 in (5.1) can be directly bounded by using the Kirchhoff 
condition (2.12) with the definitions (4.1)-(4.3), 

h-^^\\Vwh - Mo,K = h-^^\\V{w - wh) -{13- /3J||2_^ 

<\\\{w-Wh,f3-f3h)\\\h,K. (5.57) 

Step 2. We next bound the terms of rf^ in (5.2). Given now E S I{K), an edge 
of the the element K lying in the interior of 17, let 

iPE = UE{lm{(3,,)nl) , (5.58) 
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where, with a little abuse of notation, the operator 11^ is intended as appHed on each 
single component. Then, from (5.47) with integration by parts, it follows that 

= Ch]i\{LI3h,'PE)o.^ + WJ, VyjBU) , (5.59) 

where we recall that u)e was defined at the start of this section. Integration by parts 
and the equation (2.14) immediately lead to the identity 

(m(/3), V(Pe)oje = -{Lp, ^e)'^e = (9, 'Pe)u,e . (5-60) 
which, appHed to (5.59), gives 

hT\\^rr^{f3HM\\lE 

< Ch'J^ ((L/3, + q, ifEh, + {m{(3,J - m(/3), V<^s)„,) 

= Ch]^^{{Lf3f, + Qf,, ipE)oJE + (9 - 9ft, 'fiE)ujE 

+ {m{f3,,)-m{l3),y^Eh,). (5.61) 

Next, we bound the three terms on the right hand side of (5.61). For the first 
term, the identity (5.7), the Cauchy-Schwartz inequality, the definition (5.58) and the 
bound (5.48) give 



1 /2 

KC 

<C\\\{w-Wh.(3-l3^,) II [m(/3 Jn] ||o,B • (5.62) 

For the second term on the right hand side of (5.61), we note that € V* and 
ip^ = in CI\(jJe- Therefore, the duality inequality and the definition (5.58) combined 
with the bound (5.48) give 

h]l'^{Q-Qh^fEhE < ^E^I|9-9hl|-l,*,a;BllV£;||l,c^E 

< C\\q - 9j|-i,*,c.E||[m(/3Jn]||o,B • (5.63) 

For the third term of (5.61), the Cauchy-Schwartz inequality, then the inverse in- 
equality and finally (5.58) combined with the bound (5.48) lead to the estimate 

h^i\m{Ph)-m{l3),V^E)o.E < C||/3 - /3J|i,,, ^'ll<^i=;l|o,a,. 

< CII/3 - l3hh,^E\\lm{l3M\\o,E . (5.64) 

Now, by combining (5.62), (5.63) and (5.64) with (5.61) it follows 

hT\\lrn{^h)n}\WE < C{\\\{w -Wh,/3- /3J|||„,^., + ||q - gjl-i,*,..^) • (5.65) 
The remaining term of rj^ is bounded with similar arguments; with the notation 

ipE=IlE{lqh-nj), (5.66) 

the identity 

- (div q, </9e)we = {q, '^Ve)u;e (5-67) 
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with (5.54) implies 

hT\\hl-r4\\lE<Ch^i^{h-nlVE)E 

< Ch^i" ((/ - A, <p£)„, + {qy^ - q, V^e)u.^) . (5.68) 

Finally, we note that V<^b S V* and V<^b = in f2\a;£;. Therefore, 

4^119/. • rxlllcE < C(||q - + h\\\f - h\W..) ■ (5.69) 

Step 3. Third, we bound the only term of rjg e (5-3) which appears in ry^ ^ 
well. Given now a triangulation edge E in S{K) U F{K), let 

= nij(m„„(/3;,)). (5.70) 

Due to (5.47) and (2.19), integration by parts gives (here V denotes the tensor valued 
gradient applied to a vector valued function) 

= ft.^/^(Tn„(/3,, - /3), ipEn)E 

= h^\{m{p^ - /3), V(^£n))„, + {L{(3^ - /3), <fEn)^,) , (5.71) 

where n is, as usual, the chosen normal unit vector to E. For the first term, using the 
Cauchy-Schwartz inequality, then the inverse inequality and finally the bound (5.48) 
we easily get 

h]^^{m{(3f, - /3), V{ipEn))oj^ < h]^^\\(3 - (3^\\i,^^ \\V{ipEn)\\o^^j, 

< C\\f3 - /3J|i,u..||"i„„(;3J||o,E . (5.72) 

For the second term in (5.71), recalling (2.14) we have 

WE 

= h]l'^{L(3f^ + g^, </9sn)a,E + h]l.'^{q - g^, ipEn)^^ . (5.73) 

Observing now that G V* and <fEn = in Q\ue, the two terms on the right 

hand side of (5.73) can be bounded with the same arguments used above, respectively, 
in (5.62) and (5.63). Omitting the details, we therefore get 

hE^iHf^h - l3)^^En)^^ < C{\\\{w - Wh,l3 - f3h)\\\h,u,E 

+ h - qJI-i,*,we) \\m„n{Ph)\\o,E . (5.74) 

From (5.71), (5.72) and (5.74) we get 

4^'||m„„(/3J||o,B < C{\\\iw-Wh,f3-^f,)\\\h,o.E + ||q-gj|_i,.,„,) . (5.75) 

Step 4- Finally, we bound the last term ofrjp e in (5.4). Given now a triangulation 
edge E in F{K), let 

d 

= nii;(— m„s(/3^) - • n) . (5.76) 
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Using (5.47) and recalling (2.17), we obtain 

d 

- ''E^ig-^^nsi^h - /3), Ve)e + ([q - Qh] ' n, ve)e) ■ (5.77) 
For the first term, integration by parts on the edge and simple algebra give 

= h]i\{m{p - I3h)n, V^e)e - {mnn{^ - /3J, V<pe • n)^) . (5.78) 
Using again integration by parts, the first term in (5.78) can be written as 

= 4/'(i(/3 - /3J, V^b)^, + {m{f3 - /3J, VVipE^,) ■ (5.79) 

The second term in (5.77), again due to integration by parts and recalling (2.13), is 
instead equivalent to 

hE^ilQ-Qh] ■n,'~PE)E = hy'^{q- qy^,VifE)uE 

-{fh + divq,^, lpe)cje - if - fh, ^e)^e) ■ (5-80) 

For the first term, due to (2.14) and (3.11), we now have 

/iE^(9-9;i,V<^B)a,E 

= h]i\L{fit,-l3),VipE)^^-^-{^Wh-Ph,^^E)u:^), (5.81) 



where h^,^ is the size of the triangle uje- Combining all the identities from (5.77) to 
(5.81) it follows that 

< 4^'((m(/3 - /3J, VVipE)u,E - (mnniP - Hh). V<Pb • n)^ 
--^^{^Wh- (5h,VipE)wE - {fh + divqh,ipE)o,E 

-if - fh,^E)u:E) . (5.82) 

For the second term on the right hand side of (5.82), recaUing (2.17), using the 
Cauchy-Schwartz inequality and the bound (5.75) we have 

hE^{mnn{l3 - 0h),yiPE ■ n)E < h]^'^\\mnn{/^h)\\o,E hE\\VipE\\o,E 

< C{\\\{w -Wh,f3- /9J||U,c.« + IIq - qh\\-i,*,^E)hE\\^^E\\o,E , (5.83) 
which, using the inverse inequality and the bound (5.48), gives 

h^^ininnifS - f3h),^'PE ■n)E < C{\\\{w - Wh, (3 - (3h)\\\h,ujE 
d 

+ - <lh\\-h*,WE)\\-g^mns{l3h) - Qh ■ n\\o,E ■ (5.84) 
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The remaining terms on the right hand side of (5.82) can all be bounded using the 
Cauchy-Schwartz inequality, the inverse inequality and the bounds (5.56), (5.48) as 
already shown for the similar previous cases. Without showing all the details, we 
finally get 

< C{\\\{w - Wh, (3 - l3f,)\\\h,ojE + \\Q-<lh\\-h*,>^E 
+h%\\f - M|o,if) ||^m„«(/3J - Qf, ■ n\\o,E , (5.85) 

or, trivially, 

< C(|||(w; -Wh,l3- f3^)\\\h,u,^ + \\q - gj|-i,*,a,« + h%\\f - A||o,k) • (5.86) 

Recalhng now the definitions for r]K in (5.1) and the local negative norm in (5.49), 
the proposition is proved. □ 
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